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ABSTRACT 

We present here the results of our observations of TrES-2 using the Infrared Array Camera on 
Spitzer. We monitored this transiting system during two secondary eclipses, when the planetary 
emission is blocked by the star. The resulting decrease in flux is 0.127%±0.021%, 0.230%±0.024%, 
0.199%±0.054%, and 0.359%±0.060%, at 3.6 /mi, 4.5 /on, 5.8 /mi, and 8.0 /mi, respectively. We show 
that three of these flux contrasts are well fit by a black body spectrum with T e g = 1500 K, as well as 
by a more detailed model spectrum of a planetary atmosphere. The observed planet-to-star flux ratios 
in all four IRAC channels can be explained by models with and without a thermal inversion in the 
atmosphere of TrES-2, although with different atmospheric chemistry. Based on the assumption of 
thermochemical equilibrium, the chemical composition of the inversion model seems more plausible, 
making it a more favorable scenario. TrES-2 also falls in the category of highly irradiated planets 
which have been theoretically predicted to exhibit thermal inversions. However, more observations 
at infrared and visible wavelengths would be needed to confirm a thermal inversion in this system. 
Furthermore, we find that the times of the secondary eclipses are consistent with previously published 
times of transit and the expectation from a circular orbit. This implies that TrES-2 most likely has a 
circular orbit, and thus does not obtain additional thermal energy from tidal dissipation of a non-zero 
orbital eccentricity, a proposed explanation for the large radius of this planet. 

Subject headings: eclipses — infrared: stars — stars: individual (CSC 03549-02811) — techniques: 
photometric 



1. INTRODUCTION 

There has been a recent dramatic increase in the num- 
ber of extrasolar planets within 300 pc whose structures 
and atmospheric compo sitions can be probe d using the 
Spitzer Space Telescope (| Werner et al.ll200"4T ). These are 
the nearby transiting exoplanets. By measuring with 
Spitzer the decrease in light as one of these planets passes 
behind its star in an event known as a secondary eclipse, 
we can estimate the flux emitted by the planet. De- 
tection of this emission from planetary atmospheres is 
made possible by taking advantage of the enhanced con- 
trast between stars and their planets in the infrared 
wavelengths observable with Spitzer. Combining sev- 
eral of these flux measurements allows us to character- 
ize the shape of the planet's emission spectrum, which 
tells us about the properties of its dayside atmosphere. 
(Fo r a discussion of extraso l ar planetary atmo spheres, 
see lCharbonneau et all 120071 : iMarlev et al.ll2007l ). 

The planets HD 209458b and HD 189733b have been 
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the optimal choices for Spitzer studies of extrasolar 
planetary atmospheres, because of their early discov- 
ery and the relative brightness of their stellar hosts. 
We have detected infrared e mission from these exo- 
planets, both photom e trical ly (|Deming et al.l[2005l . 120061 : 
iKnutson et all 120071. l2008h. and wit h low-resolution 
spect ro scopy ( Richard son et alj 120071: iGrillmair et alj 
2007: iSwain et alj |200? ~ 
Grillmair^taLl 12008 



Charbonneau et al.l 2008; 

An- 



_, ISwain et al.l 120081. 12009, 
other early success was the report bv lCharbonneau et al.l 
(2005) of the m easurement of infr a red lig h t from TrES-1. 
More r ec ently. Harrington et al.l (120071 ) . iDeming et al 
(|2007bf) . iDemorv et alj (|2007t ). and iMachalek et~a! 
(2008) announced the results from Spitzer observations 
of the transiting exoplanets HD 149026b, GJ 436b, and 
XO-lb. 

There has been a flurry of activity to reconcile atmo- 
spheric models with this limited number of infrared mea- 
surements. While several attempts have been made to 
expla i n the infrared observa t ions (s e e, e.g. , Barman et al 
2005; IBurrows et all 120051 . I2007al : IFortnev et al.l 1200 
120061 : ISeager et al.ll2005l ). the models are not entirely in 
agreement and n o single model can expl a in every obser 
vation. Rece ntly. IBurrows et al.l (l2007bf ). IB 
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(2008), and IFortnev et ail (|2008f ) supplied a possible 
piece of the puzzle by proposing that the very highly ir- 
radiated hot J upiters such as HP 2 09458b and TrES-2 
(see Fig. 1 of IFortnev et al.l I2008 P) will exhibit water 
emission rather than the expected water absorption at 
the IRAC wavelengths, a result of a temperature in- 
version in their atmospheres. The large stellar insola- 
tion that TrES-2 experiences may permit the presence 
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of TiO and/or VO molecules in the hot planetary at- 
mosphere that would condense in a cooler atmosphere. 
These opaque molecules would then cause the inver- 
sion. However, while this hypothesis can explain the 
emission features seen in the spectrum o f HD 209458b 
(jBurrows et al.ll2007bl ; iKnutson et al.ll2008l ). it is incom- 
plete. The relatively lower insolation experienced by the 
planet XO-lb is similar to that of TrES-1, and hence both 
atmospheres would be predicted to have water absorp- 
tion features, as is inde ed supported by the infra red ob- 
servations jjDfJTrJ^lbyJCh^ ([2005) . Nev- 
ertheless, Machalck et al. (2008) have shown that XO-lb 
displays the contrary, with evid ence of water e missio n in 
its atmosphere. More recently, iFressin et al.1 ()2009f ) did 
not find any evidence for the expected thermal inversion 
in the atmosphere o f the highly irradiated planet TrES-3. 
IZahnle et al.l (|2009T ) have shown that sulfur photochem- 
istry may also lead to the formation of inversions. In 
their models, this photochemistry is more or less tem- 
perature independent between 1200 K and 2000 K, and 
is also relatively insensitive to atmospheric metalicity. 

Although over 400 extrasolar planets are known, it is 
only for these nearby transiting exoplanets that we can 
measure the planetary radii and true planetary masses 
precisely enough to provide useful constraints for the- 
oretical models. There have been problems reconcil- 
ing the observed planetary masses and radii with mod- 
els (see lLaughlin et al.l 120051. iCharbonneau et al.l 120071 . 
iLiu et al.ll2008l . and llbgui & Burrowsll2009l . for areview) , 
namely that there arc planets such as TrES-2 whose 
radii are larger than predicted by the standard mod- 
els for an irradiated gas giant planet. Several explana- 
tions have been proposed for bloated planets like TrES-2, 
mainly related to some additional source of energy to 
combat planetary contraction. One such possible energy 
source is the tidal damp i ng of a non-zero eccentricity 
(jBodenheimer et al.l l200ll . 120031 ): the planetary orbit of 
a hot Jupiter is expected to be circular, unless it is grav- 
itationally affected by an uns een planetary companion 
(see, e.g., iRasio fc Fordl 119961 ) . We can constrain the 
likelihood of a tidal damping energy source by measur- 
ing the timing of the secondary eclipse of a transiting hot 
Jupiter like TrES-2 and comparing these timings to pre- 
dictions based on transit timings a nd the hypothesis o f 
a circular planetary orbit (see, e.g.. lDeming et aLlfe OOS). 

Th e transiting hot Jupiter TrES-2 ( |Q 'Donovan et all 
2006) is one of the known hot Jupiters with a radius 
larger than expected from current models. The atmo- 
sphere of TrES-2 experiences similar levels of irradia- 
tion from its host star as the atmosphere of the likewise 
bloated planet HD 209458b, and hence we expect to find 
evidence of a ther mal inversion for TrES - 2, according to 
the pr edicti ons of IBurrows et al.l (|2007bl ) , iBurrows et al.1 
( 2008), and lFortnev et al.l (|200l ). Here we present the 
first Spitzer observations of TrES-2 (*J1]). From our analy- 
sis (SJ3|, we have detected thermal emission from the tran- 
siting planet. We found no evidence for timing offsets of 
the secondary eclipses, and deduce the possible presence 
of a thermal inversion in the atmosphere of TrES-2 (fQJ. 

2. IRAC OBSERVATIONS OF TRES-2 

We monitored TrES-2 using Spitzer during the 
time of two secondary eclipses, employing a differ- 
ent pair of the four bandpasses available on the In- 



TABLE 1 
TrES-2 System Parameters 



Parameter 



Value 



Reference 



P (d) 
T c (HJD) 
b = a cos i/Ft* 
(°) 

Rp/R* 

T s [3.6 Atm/5. 8^111] (HJD) 
T B [4.5Atm/8.0Atm] (HJD) 



2.470621 ±0.000017 
2,453,957.63479 ± 0.00038 
0.8540 ± 0.0062 

83.57 ± 0.14 
0.1253 ±0.0010 

2,454,324.5220 ± 0.0026 
2,454,070.04805 ± 0.00086 



a IHolman etTaTl 1120071) . 
b This work. 



frared Array Camera (IRAC; iFazio et all [2004) dur- 
ing each of the eclipses. We took care to position 
TrES-2 (2MASS J19071403+4918590: J = 10.232 mag, 
J — K s = 0.386 mag) away from array regions impaired 
by bad pixels or scattered light. We also kept the cor- 
responding IRAC stray light avoidance zones free of 
stars that are bright in the infrared. We observed a 
5f2 x W.2 field of view (FOV) containing TrES-2 dur- 
ing two eclipses, using the Stellar Mode of the IRAC 
instrument. On UT 2006 November 30 (starting at 
HJD 2,454,069.956), we obtained 1073 images of this 
FOV at 4.5 /jm and 8.0 fim. with an effective integration 
time of 10.4 s for a total observing time of 3.9 hr. Our 
3.6-/im and 5.8-/zm observations of TrES-2 were taken 
on UT 2007 August 16. The observations began at 
HJD 2,454,324.436 and lasted 4.0 hr, during which we ac- 
quired 2130 and 1065 images in the respective channels 
(of effective exposure time 1.2s and 10.4s, respectively). 

3. DERIVING AND MODELING LIGHT CURVES OF TRES-2 

As part of the Spitzer Science Center pipeline for IRAC 
data (version S15.0.5 for the 2006 observations, version 
S16.1.0 for the 2007 observations), the images were cor- 
rected for dark current, flat-field variations, and certain 
detector nonlinearities. Each header of these Basic Cali- 
brated Data (BCD) images contains the time and date of 
observation and the effective integration time. We used 
these to compute the Julian date corresponding to mid- 
exposure of each observation. In order to convert these 
dates to Heliocentric Julian dates, we calculated the cor- 
responding light travel time between the Sun and Spitzer 
(ranging from 0.5 to 2.5 minutes), using NASA JPL's 
HORIZONS^ service. We then computed the orbital 
phases using an updated orbital period an d transit epoch 
(see Table [T]) derived from observations ( Holman et al. 
120071 ) made as part of the Transit Light Curve project. 

Using the nearest integer pixel (xoj2/o) as an initial 
estimate of the position of TrES-2 on the array, we com- 
puted the flux- weighted centroid (x, y) of TrES-2 in each 
BCD image. The intra-pixel position of TrES-2 is then 
(x' = x — xq, y' = y — yo). We measured the flux from 
our target using circular apertures ranging from 2 to 10 
pixels, and subtracted the background signal assessed in 
a sky annulus with inner and outer radii of 20 and 30 
pixels, respectively. We normalized the fluxes for a given 
channel and aperture size by dividing each time series 
by the median of its values outside the times of eclipse. 

8 http : //ssd . jpl .nasa. gov/?horizons 
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We examined the variation of the residual scatter of the 
out-of-eclipse data, and found that an aperture radius 
of 3 pixels produced the smallest residual scatter for all 
four channels. Figure Q] shows the four light curves we 
obtained using this photometric aperture. 

For each of the observed secondary eclipses, we cre- 
ated a model of the two IRAC light curves obtained dur- 
ing that eclipse. We first accounted for various detector 
effects known to be present in IRAC data. There is a 
known correlation between the IRAC 3.6-/im or 4.5-/Ltm 
flux from a sour ce and the intra-p i xel position (x',y') 
on the detector (iReach et all 120051 : iCharbonneau et all 
120051 . 120081 : Knutson et all 120081 ): the sensitivity of an 
individual pixel varies depending on the location of the 
stellar point-spread function (PSF), with higher fluxes 
measured near the center of the pixel and lower fluxes 
near the edges. Figure [1] shows this effect in our TrES-2 
data. Data from these two channels also demonstrate 
a linear trend with time (as previously observed by 
iKnutson et ail 120091 in IRAC observations of TrES-4), 
with a positive trend at 3.6 (im and a negative slope at 
4.5 /im. We removed these trends to obtain /o, the actual 
stellar flux: 

fo = f/([ci + c 2 x' + c 3 y' + c 4 (x') 2 + c 5 (y') 2 } x [1 + Cdt}), 

where / is the measured stellar flux, (x' , y 1 ) is the intral- 
pixel position, dt is the amount of time from the first 
observation, and (ci_5,C) are free parameters in our 
model. The 5.8-psm and 8.0-/zm data (see Fig. [1]) display 
the "ramp" associated with th es e IRA C channels first 
noticed bv ICharbonneau et all (|2005l l and expanded 
upon bv lHarrington et al.l (|2007l ) (supplementary infor- 
mation). Both data sets showed an overall increase in 
flux with dt. For these two data sets, we removed this 
detector effect by including the following correction in 
the model: 

fo = f/[dx + d 2 (\ndt') + d 3 (\ndt') 2 ], (2) 

where d±- 3 are free parameters, and dt' — dt + 0.02 days. 
Here our substitution of dt' for dt is simply to avoid divi- 
sion by infinity. We then modeled simultaneously the two 
corrected time series using the eclipse light curve equa - 
tion for a uniform source from IMandel fc Ago] ([2002). 
We obtaine d the required syste m parameters (see Ta- 
ble [U from iHolman et ail ([2007): the planetary orbital 
period, impact parameter, orbital inclination, and the ra- 
dius ratio between the planet and the star. Based on this 
ephemeris, we calculated the predicted eclipse epoch (T s ; 
see Table [1]). The three free parameters for the eclipse 
model were the timing offset (At), the eclipse depth at 
the shorter wavelength (A/;), and the depth at the longer 
wavelength (Afh). The model of the two light curves 
therefore has 12 free parameters: [c\, C2, c 3 , C4, C5, C, 
di, d 2 , d 3 , At, Af u Af h ]. 

For a particular instance of this model, we compute 
the x 2 measure of the goodness of fit to the two relevant 
data sets as follows. For each light curve, we compute an 
initial uncertainty for the normalized flux values as their 
standard deviation. We exclude 5-ct outliers in flux from 
further consideration. We then exclude large outliers 
in intra-pixel position (max [\x' — x' m \, \y' — y' m \] > 0.15, 
where [x' m ,y' m ] are the median intra-pixel values) on the 
detector. We compute the \ 2 01 the two light curves 



TABLE 2 
Best-Fit Values for Free 
Parameters of Eclipse Models 



Parameter 


Value 


3.6 /.irn/5.8 /.im 




x 2 


3164 


N 


3166 


At (min) 


1.8 ±3.6 


A/,( = A/3.6 „m) 


0.127%±0.021% 


Afh(= A/5.8 fim) 


0.199%±0.054% 


4.5 fim/S.O /im 




x 2 


2117 


N 


2119 


At (min) 


0.7±3.1 


A/;(= A/4.5 „m) 


0.230%±0.024% 


Af h {= A/g.O^m) 


0.359%±0.060% 



separately, and rescale each \ 2 to reflect the reduction of 
the number of data points by the above outlier exclusion. 
We then sum the resulting values to compute the \ 2 of 
the overall model. 

To find an initial estimate of the best-fit parameters 
for each model of t wo light curves, w e first used the 
AMOEBA algorithm (jPress et al.lll992l ) to minimize the 
X 2 of the fit. Using this initial estimate as a starting 
point, we applied the Markov Chain Mon te Carlo method 
(see, e.g.. iFordl 120051 : iWinn et all 120071 ). computing the 
X 2 at each of the 10 6 steps of the chain. We then calcu- 
lated the median of the 10 6 x 2 values, and excluded from 
further analysis all the steps prior to the occurrence of 
the first value lower than this median. For the ith free 
parameter of our model, we derived the best-fit value pi 
as the median of the remaining distribution of values for 
that parameter. We computed the (possibly unequal) 
lower (er-,i) and upper (cr+,i) errors in this value such 
that the ranges [pi —a^ } i,pi] and [pi,Pi + cr+ t i] each con- 
tain 68%/2 of the values less than or greater than, respec- 
tively, the best-fit value. We computed the flux residuals 
after dividing out the best-fit model, and then derived 
an updated uncertainty as the new standard deviation of 
the residuals from the model. We again computed and 
rescaled the x 2 °f the two light curves separately, and 
then summed the resulting values to compute the x 2 of 
the best-fit. 

We tabulate for both pairs of channels the best-fit val- 
ues for the free parameters of the eclipse models and the 
reduced x 2 m Table O We have overplotted in Figure Q] 
the above corrective functions using the best-fit parame- 
ters we derived from the Markov chains. In Figure^ we 
plot the corrected fluxes from TrES-2 at the four wave- 
lengths, and overplot the best-fit eclipse models. In both 
figures, the error bar shown for each binned data point 
is the standard deviation of the flux values in that bin, 
divided by the square root of the number of points in the 
bin. 

4. DISCUSSION AND CONCLUSIONS 

From the best-fit values (see Table [5]) for the timing 
offset of the two observed secondary eclipses, we see that 
their weighted average (At = 1.2 ± 2.3 min) is consistent 
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Fig. 1. — Relative fluxes from the system TrES-2 at 3.6 fim, 4.5 fim, 5.8 /im, and 8.0 fim (with an arbitrary flux offset), binned and plotted 
vs. the time from the predicted center of secondary eclipse (T s ). Superimposed are our best-fit models (black lines) for known instrumental 
effects. The error bar shown for each binned data point is the standard deviation of the flux values in that bin, divided by the square root 
of the number of points in the bin. 



with no offset from the predicted epochs for the eclipses. 
Note that we have no t accounted for the light travel delay 
time fsee lLoebl [2005( 1 of 37 s across the TrES-2 system, 
because of the relatively large size of the errors in these 
timing offsets compared to this delay time. An upper 
limit for the orbital eccentricity of a transiting planet can 
be computed from the timing offset At, using e cos a; ~ 
7r At/2 P, where uj is the unknown lo ngitude of periastron 
and P is the known orbital period (|Charbonneau et al.l 
120051 . Eq. 4). The 3-cr upper limit for TrES-2 is therefore 



0.0036, consistent with a negligible orbital eccentricity, 
unless w ~ 90°. Tidal dampin g of orbital eccentricity 
(|Bodenheimer et al.ll200i1 . 12003D is therefore unlikely to 
be a sizable contribution to the internal energy of this 
bloated exoplanet. 

We now turn to a discussion of the TrES-2 planet- 
star contrasts. We first emphasize how well the data 
can be fit by a blackbody spectrum with no molecular 
band features. The black filled circles with error bars 
in Figure [3] show the Spitzer IRAC data points from 
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Fig. 2. — Same relative fluxes as Figure^ except that here the fluxes have been corrected for the known detector effects. Overplottcd 
are our best-fit models (black lines) of the secondary eclipses. The error bars are as in Figurefl] 



this work. The blue dashed line is a 1500 K blackbody 
flux divided by the Kurucz stellar model with stellar pa- 
rameters (T cff = 5750 K, Z = 0.0, log(g) = 4.5) clos- 
est to those (T ea = 5850 ± 50 K, Z = -0.15 ± 0.10, 
log(» = 4.4 ± 0.1) derived bv lTorres et~aTl ([2008T ). It is 
clear that the blackbody fits all the data points resonably 
well, given the large error bars. However, a blackbody 
spectrum is only a nominal guideline, since the actual 
planetary spectrum is influenced by the myriad contri- 
butions due to molecular band features, collision induced 
opacities, temperature gradients. 



Model atmospheric spectra for TrES-2 are also shown 
in Figure El The red and green spectra show models 
with and without a thermal inversion on the planet day- 
side, respectively. The red and green circles (enclosed in 
black circles) show the corresponding model points ob- 
tained by integrating the spectra over the Spitzer IRAC 
bandpasses. The corresponding model thermal profiles 
are shown in Figure |U The spectra were generated 
using the hot Jupiter atmosp here model developed in 
Madh usudhan fc Seageri (|2009f ). We consider a cloudless 
atmosphere, and the molecular species are assumed to be 
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Fig. 3. — Contrast ratios (black filled circles) for TrES-2 at 3.6 £tm, 5.8 /im and 8.0 fim, which are consistent with a model (blue line) 
of a 1500 K black body planetary flux divided by a Kurucz model of the star TrES-2. Our additional observation at 4.5 /im shows some 
evidence of excess emission at this wavelength. Also shown are the predictions (red and green ci rcles) for these four fluxes from theoretical 
planet-star flux contrast models (red and green lines) computed for the star TrES-2 using the Madhusudhan & Scagcr (2009) code (see 
text). The red (green) spectrum corresponds to a model with (without) a thermal inversion on the planet dayside. 



well mixed. The stellar spectrum was represented by the 
appropriate Kurucz model. The model spectrum with- 
out a thermal inversion has T e ff = 1634 K, and the model 
spectrum with a thermal inversion has T e g = 1459 K. 
Both the models allow for extremely efficient day-night 
energy redistribution. At face value, we find that both 
the models fit the data almost equally well. 

However, the two models require different molecular 
compositions, which helps us determine the more prob- 
able model. The model without a thermal inversion has 
uniform molecular mixing ratios of 10~ 4 for H2O, 1CP 6 
for CO, and 10 -6 for CH4. While the mixing ratio of 
H 2 is plausible, the low mixing ratio of CO is surpris- 
ing. In a hot atmosphere, with T c h = 1500 K, CO is 
expected to be highly abundant. On the other hand, the 
model with a thermal inversion has uniform molecular 
mixing ratios of KT 4 for H 2 0, 10~ 4 for CO, 5 x 10~ 5 for 
CH4, and 2 x 10 -6 for CO2. These compositions show 
a relatively high abundance of CO, as expected in a hot 
Jupiter atmosphere like that of TrES-2. Based on these 
compositions, the presence of a thermal inversion seems 
like a more favorable scenario. Additionally, the require- 
ment of 2 x 10~ 6 of CO2 for the inversion m odel suggests 
an enhanced metallicity of [M/H] ~ 0.7 (jZahnle et al.l 



2009). However, a thorough exploration of the parame- 
ter space is needed to place constraints on the thermal 
inversion in conjunction with th e molecular compositions 
(jMadhusudhan fc Seagerll2010l ). 

Thus we see some evidence that TrES-2 is part of the 
group of planets that display excess emission at wave- 
lengths longer than 4 /im. The presence of a thermal in- 
version i n the atmosph e re of TrE S-2 would support the 
theory of lFortnev etail (j2008f ) and lBurrows et al.l (|2008l ) 
that planets with substellar fluxes greater than or equal 
to approximately 10 9 erg s _1 cm~ 2 should show evidence 
for a temperature inversion in their atmospheres. 

Our theoretical models for the planetary flux from 
TrES-2 are currently constrained only at infrared wave- 
lengths. Such an infrared s ignal appears to be r elatively 
insensitive to cloud cover ([Burrows et al.ll2005l ). How- 
ever, the amount of starlight reflected off the planetary 
atmosphere (at optical wavelengths) is highly dependent 
on the presence and size of upper-a tmospheric conden- 
sates such as MgS iOs, Fe and A1 2 3 (jMarlev et al.lll999l ; 
iGreen et al.ll2003t ). Hence, a tighter constraint of atmo- 
spheric models for TrES-2 would be derived from the 
combination of measurements of the optical planetary 
flux from observations of secondary eclipses with the re- 
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suits from our infrared observations. Two likely sources 
for such observations of the optical flux fro m TrES-2 
are the Kepler mission (iBorucki et al.l I1997T ) and the 
EPOXI mission (jDeming et al.l l2007af h Kepler has al- 
ready demonstrated the ability to monitor the variation 
in thermal emission and reflected light from HAT-P-7 
(|Borucki et al.ll2009l ). We eagerly await the comparison 
of our infrared observations with future optical observa- 
tions from these two missions. 

Although our observations of the atmospheric emis- 
sion from TrES-2 at 3.6 /im, 4.5 /im, 5.8 ^m, and 8.0 /im 
provide limited spectral coverage, we have been able to 
deduce the probable presence of gaseous molecules with 
high opacities in the atmospheric that result in emis- 
sion in the 4.5 ^m band. The highly irradiated gas giant 
TrES-2 thus may provide additional evidence for the cor- 
relation between the occurrence of thermal inversion in 
the atmosphere of a planet and the level of stellar inso- 
lation experienced by the planet. 

Our observations of two secondary eclipses by this ex- 



oplanet occurred at the time predicted using the time of 
transit and the assumption of a circular orbit, within the 
errors. From this we conclude that the source of addi- 
tional energy required to inflate the planetary radius to 
its bloated size is unlikely to be tidal heating caused by 
the circularization of an eccentric orbit. 



This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under a contract with the National Aeronautics and 
Space Administration (NASA). This research was sup- 
ported in part by NASA under grant NNG05GJ29G (is- 
sued through the Origins of Solar Systems Program) , and 
also by an appointment to the NASA Postdoctoral Pro- 
gram at the Goddard Space Flight Center (administered 
by Oak Ridge Associated Universities through a contract 
with NASA). 
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Fig. 4. — Pressure- Temperature profiles and contribution functions corresponding to the models shown in Figure [3] The left panel 
shows the two pressure temperature profiles (with and without a thermal inversion) corresponding to the two models reported in this 
work. The solid (dash-dot) curve shows the P-T profile with (without) a thermal inversion. Both the profiles lead to good fits to the 
o bservations, albeit with d i fferen t molecular compositions. The P-T profiles were generated with the parametric prescription developed 
in Madhusudhan & Scagcr (2009). The thermal inversion in the solid profile lies between 0.20-0.01 bar, and spans temperatures between 
1100-2000 K. The right panel shows the contribution functions in the six Spitzer photometric channels for each model. The legend shows 
the channel center wavelength in microns, and the curves are color-coded by the channel. The solid curves show the contribution functions 
for the inversion model, whereas the dash-dot curves show the contribution functions for the non-inversion model. All the contribution 
functions are normalized to unity. 



